Viscosity of strongly interacting quantum fluids: spectral functions and sum rules 
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The viscosity of strongly interacting systems is a topic of great interest in diverse fields. We focus 
here on the bulk and shear viscosities of non-relativistic quantum fluids, with particular emphasis on 
strongly interacting ultracold Fermi gases. We use Kubo formulas for the bulk and shear viscosity 
spectral functions, ({uj) and r?(tj) respectively, to derive exact, non-perturbative results. Our results 
include: a microscopic connection between the shear viscosity rj and the normal fluid density p„; 
sum rules for ({ijj) and r]{Lj) and their evolution through the BCS-BEC crossover; universal high- 
frequency tails for and the dynamic structure factor ^(q, oj). We use our sum rules to show 
that, at unitarity, ({ijj) is identically zero and thus relate r;(tj) to density-density correlations. We 
predict that frequency-dependent shear viscosity r;(tj) of the unitary Fermi gas can be experimentally 
measured using Bragg spectroscopy. 

PACS numbers: 67.85. De, 67.10.Jn, 67.85. Lm 



C/3 



I 

o 



(N 
> 
OS 

00 

o 

O 

o 



X 



I. INTRODUCTION 

The study of the viscosity of strongly interacting quan- 
tum fluids has brought together very different areas of 
physics - black holes and string theory, quark-gluon plas- 
mas, quantum fluids and cold atoms - which, at first 
sight, appear to have little in common ^2,]. This ex- 
traordinary development originated with the work of Son, 
Starinets and coworkers (T], y, 01 who calculated the shear 
viscosity in a strongly interacting quantum field theory, 
the A/" = 4 supersymmetric Yang-Mills (SYM) theory, 
and conjectured a lower bound 



i]/s > h/{'i-KkB) 



(1) 



for the ratio of the shear viscosity r] to the entropy den- 
sity s of any system. These results were obtained using 
the AdS / CFT formalism where certain strongly coupled 
field theories can be mapped onto weakly coupled gravity 
theories. 

Although a number of counterexamples have since 
been proposed 0-[l|> there are no known experimen- 
tal violations of the bound given by Eq. ([1]). Remark- 
ably, two very different experimental systems come close 
to saturating the bound: the quark-gluon plasma at 
Brookhaven's Relativistic Heavy Ion Collider J), 10], and 
ultracold atomic Fermi gases [ll|, [l^j close to a Feshbach 
scattering resonance, where the s-wave scattering length 
becomes infinite 13[. This is the strongly interacting 



unitary regime that lies at the center of the BCS-BEC 
crossover. These two systems are amongst the hottest 
and coldest systems every realized in a laboratory. 

In this paper we focus on non-relativistic quantum 
fiuids, with particular emphasis on strongly interacting 
Fermi gases. These are systems for which the most 
controlled experiments should be possible. A "perfect 
fluid" with the minimum shear viscosity is necessarily in 
a quantum regime, since the bound involves h. In ad- 
dition, it must also be in a strongly interacting regime 
where well-defined quasiparticle excitations do not exist. 



If the system had sharp quasiparticles, then their mean 
scattering rate r^^ would be much less than the aver- 
age energy per particle eo, so that h/r <^ e^. We can 
then use Boltzmann's kinetic theory approach to obtain 
rj ~ neoT, where n is the number density. Using s ^ nks, 
we flnd a large r]/s ^ t^T/kB ^ h/ks. Thus, in order 
to find perfect fluids that come close to saturating the 
lower bound given by Eq. ([1]), one must look at strongly 
interacting quantum fluids where the quasiparticle ap- 
proximation fails. 

In this paper we use Kubo formulas for the frequency- 
dependent spectral functions for shear viscosity ?y(a;) and 
bulk (or second) viscosity C('^)j and derive several exact, 
non-perturbative results without making weak coupling 
or quasiparticle approximations. Our main results are: 

• We establish a microscopic connection between the 
shear viscosity r] and the normal fluid density p„ and 
show that a non-zero p„ is a necessary condition for a 
non- vanishing r]. 

• We derive sum rules for rj{u)) and Q{uj) of any Bose 
or Fermi system with an arbitrary two-body interaction; 
see Eqs. ^ and ((55)) . 

• For a dilute two-component Fermi gas, we find the 
shear viscosity sum rule 



duj 



C 



C 



(2) 



valid for arbitrary temperature and l/(A;i7'a), where a is 
the s-wave scattering length. Here, e is the energy den- 
sity and C is the contact [14| . A central quantity in many 
of our results, C = k%C[l/ {kFa),T / ep] can be defined via 
the large-fc tail of the momentum distribution rik ~ C/Zc^ 
for A: ^ kp, and characterizes the short-distance proper- 
ties of the many-body state. 



For the bulk viscosity, we obtain the sum rule 



dijj C,{uj) 



1 



727rma2 V da 



dC 



(3) 



2 



where the derivative is at fixed entropy per particle s = 
S/N. [Different from Eq. ([T]), in the remainder of this 
paper we use s to denote this quantity rather than the 
entropy density.] Below the superfluid transition, the 
bulk viscosity that enters Eq. ([3]) is C2, associated with 
the damping of in-phase motions of the superfluid and 
normal components. The positivity of C(<^) implies that 
idC/da-^)s > 0. 

• At unitarity, the bulk viscosity spectral function van- 
ishes at all frequencies and all temperatures. Quite gen- 
erally, C{uj) > 0, but the sum rule in Eq. ([3|) vanishes 
for \a\ — OD and thus Ci'-^) — for the unitary Fermi 
gas. This generalizes the result [HI that the static bulk 
viscosity ^(0) vanishes at unitarity. 

• It follows from the previous result that, at unitarity, 
the shear viscosity spectral function ri{Lo) can be related 
to density-density correlations as 

T]{uj) = linr — Imxpp(q,a;) (|a| = 00). (4) 

Thus, 77(0;) for the unitary Fermi gas can be measured 
spectroscopically using, for instance, two-photon Bragg 
spectroscopy. 

• We show from our sum rules that various dynamic re- 
sponse functions for Fermi gases have high-frequency tails 
characterized by odd-integer power laws, whose magni- 
tudes are controlled by the contact C. The tail C/^/uj 
of r](uj) is evident from Eq. Using this we find that 
the dynamic structure factor has a tail of the form [l^ 
5(q,cj) ~ Cg4/a;7/2 for g ^ and a; -> 00, which is 
shown to be a generic feature of short range physics. 

In the remainder of this Section we describe how the 
rest of the paper is organized. In Section [III we begin 
with a careful derivation of Kubo formulas for the spec- 
tral functions rj^tj) and C(j-^) in terms of current-current 
correlation functions, Eqs. (|20|) and ((2T|) . and, equiva- 
lently, in terms of the stress-stress correlator, Eq. ([25]) . 

In Section IIIII we recall some elementary facts about 
the shear viscosity of a fluid and why it is analogous to 
the resistivity, and not the conductivity, of a metal. We 
then establish a connection between the viscosity 77 and 
the normal fluid density pn using microscopic response 
functions. 

After establishing the positivity of ri{uj) and of ({lu) 
in Section IIVI we derive sum rules for these quantities 
in Section |Vl The most general sum rules for the shear 
and bulk viscosities of any Bose or Fermi system with an 
arbitrary isotropic interaction potential V{p), and valid 
for all temperatures, are given in Eqs. (|52p and ([5^. 

In Section |Vll we specialize to the dilute Fermi gas, 
with interparticle spacing kp^ and s-wave scattering 
length a both much larger than the characteristic range 
To of the potential. We obtain the ( sum rule in Eq. ([3]) , 
which is finite in the zero range limit. The rj sum rule, 
however, has an ultraviolet divergence; see Eq. (1551) . We 



identify, in Section [VIIl the C/^/w high-frequency tail 
of the shear viscosity spectral function, and derive the 
sum rule given by Eq. ([2|), which is manifestly finite for 
To — 0. The sum rules given by Eqs. ([2]) and ^ are 
valid in both the normal and superfluid phases, with ^ 
replaced by (2 in the latter state. 

In Section IVlIII we show from the C sum rule that, at 
unitarity, Ci'-^) vanishes at all frequencies and all temper- 
atures. We also discuss the l/(fci?a)-dependence of the 
7] and C sum rules across the BCS-BEC crossover, us- 
ing available quantum Monte Carlo data for the energy 
density at T = 0. 

We discuss the connection between viscosity and 
density-density correlations in Section IIXI and find two 
interesting results. First, we show how a density probe 
such as two-photon Bragg spectroscopy can in principle 
be used to measure the shear viscosity spectral function 
r]{uj) at unitarity. Second, we identify the high-frequency 
tail in the dynamic structure factor S{q,uj). 

In Section |Xl we briefly compare the sum rules that 
we have derived for non-relativistic quantum fluids with 
those obtained in relativistic quantum field theories. Fi- 
nally in Section [XII we conclude with open questions. 

There are five Appendices which contain technical de- 
tails of derivations or review certain results which are 
used at various places in the paper. In Appendix [XJ 
we briefly discuss an alternate stress tensor operator of- 
ten used to calculate the shear viscosity. Some results 
from dissipative two-fluid hydrodynamics, which we use 
in our paper, are reviewed in Appendix IB] We review in 
Appendix [C] results related to the contact that are used 
at several places in the paper, and also give a detailed 
derivation of certain equations that involve the contact. 
In Appendix [Dj we derive a microscopic expression for 
the pressure. Finally, in Appendix |E] we give details of 
the derivation of the C sum rule which make use of the 
scaling form of thermodynamic functions across the BCS- 
BEC crossover. 



II. KUBO FORMULA FOR VISCOSITY 

We begin by deriving Kubo formulas for the bulk and 
shear viscosity. Although the results of this Section are, 
for the most part, "well known", we could not find a 
complete derivation at any one place in the literature. 
In particular, there are several subtle points not dealt 
with adequately elsewhere, not least the definition of the 
stress tensor operator Hap for non-relativistic systems. 

To introduce notation, we start with the Euler equa- 
tion 

mdtja{r,t) ^ -di3llai3{r,t), (5) 

where m is the mass of the particles, ja is the (number) 
current and 11^^ is the momentum flux density tensor, 
which we call the stress tensor, for short. Here, a and /3 
take on values x, y, z (and there is no difference between 



3 



upper and lower indices in our non-relativistic formula- 
tion). In general, the stress tensor is given by fl7{ 

^aP = PSal3 + pUaUfj - a'afj, (6) 

where P is the pressure, p the mass density and u the 
velocity. The viscous term a'^p is given by 



dfjUa+daUp- -SapiV ■ u) 



-C<5„^(V-w), (7) 



where rj is the shear viscosity and C the bulk viscosity. 
The generalization of Eq. ([7]) to the superfluid state is 
well known and involves additional bulk viscosities. 
At the end of Section III Al we show that the Kubo for- 
mula we derive for ( describes the bulk viscosity ^2 in the 
superfluid phase. 

Our goal is to obtain Kubo formulas for frequency- 
dependent generalizations of the long- wavelength viscosi- 
ties, r] and C: in terms of equilibrium correlation func- 
tions of the many-body system. The Kubo formulas for 
viscosities are often written in terms of the stress-stress 
correlators; see, e.g.. Sec. 90 of Ref. [11]. However, the 
form of the stress tensor (or momentum flux density) 
operator is not obvious, and many different, compli- 
cated expressions [l^ which are presumably equivalent 
can be found in the literature. Part of the problem is to 
write down an operator expression for the pressure P in 
terms of particle positions and momenta. In high-energy 
physics, a simple way to calculate the stress-energy ten- 
sor Hafj is to vary the action with respect to the metric in 
curved space-time. We prefer, however, to describe non- 
relativistic fluids without going to curved space-time. 

To begin with, in III Al we adopt an approach that per- 
mits us to get around the complexities of defining the 
stress operator 11. We consider the linear response of 
a fluid to an externally imposed velocity field and derive 
Kubo formulas for the bulk and shear viscosities in terms 
of current- current correlators. The results of this subsec- 
tion are the same as those of Kadanoff and Martin f20| . 

In lllBI we use an operator form of Euler's equation to 
make the connection between bulk and shear viscosities 
and stress-stress correlators. In Appendix \^ we derive 
an alternative form of the stress correlator, which works 
only for the shear viscosity in the zero- frequency limit, 
but is often used in calculations. 



A. Current correlators 

We calculate within linear response theory [2l|, the 
current flow in a fluid subjected to an external velocity 
field u{r,t) = M(r)e~*'^*e° * which is turned on adiabat- 
ically. Our goal is to relate the imaginary part of this re- 
sponse function to viscosity through the dissipative part 
of the stress tensor. 

The response of a fluid to the "moving walls" of its 
container is a standard concept in the theory of super- 
fluidity 22|. Here, we generalize this analysis to a non- 
uniform and time- varying external perturbation u(r,t), 



taking the long wavelength limit at the end. We write the 
Hamiltonian of the system H plus external perturbation 
H' as 23] 



N 



total 



N 

= H - - Y,j dvu{v,t)-{v,,5{v~i,)} + 0[u%{8) 



where pi and f ^ are the momentum and position opera- 
tors, respectively, for the «th particle, m is the mass, and 
V is the potential energy operator. The anticommutator 
{A, B} — AI3 + BA is used to symmetrize products. We 
thus see that to linear order in u, the external perturba- 
tion is 



H'{t) = -m I dre-*"*e°''*M(r) •j(r,t), (9) 



where j = X^iLi {P*; 5{r~Ti)} /2m is the current density 
operator. 

Linear response theory gives the result [2l], 



Oa(r,i)) =m dr' 



dt'e 

— oo 



Xf{r~v\t^t')up{r'). 



(10) 



Here and below, we use the standard convention of sum- 
ming over repeated indices. The retarded current cor- 
relation function Xj^ is obtained by using A — ja and 
B = in Eq. ^ below. 

For later use, we provide a general deflnition for the 
retarded response function, or correlator, for operators A 
and B : 



XA,B(r-r',<-i') = 

iQit~t'){[A{r,t),B^iv',t')]). 



(11) 



Here, (Q) = Tt[Q exp{-H /T)]/ Z is the thermal expec- 
tation value at temperature T and Z — Trjexp(— if/T)] 
is the partition function. The step- function &{t — t') en- 
forces causality. We will use the convention of unit vol- 
ume = 1 and set h = = 1, unless explicitly stated 
otherwise. 

We flnd the spectral representation for xa,b using the 
exact eigenstates and eigenvalues of the fully interact- 
ing many-body Hamiltonian H\a) = Ea\a), and Fourier 
transform the result to obtain 

XAs(q,c^)-^5I^~^"^^^ 

a.b 

{a\Bl\b){b\A^\a) {a\A^\b){b\Bl\a) 



(12) 



UJ + Eba+ iO+ UJ-Eba+ iO+ 

where Eba = Eh — Ea. The quantity of central interest 
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to us in this paper is the imaginary part of given by 
ImxA,B(q,^) = 7r(l-e-"/^) 

X I E (^l^ql^) (b\Bl\a)6{u; - E,,). (13) 

a,b 

Returning to the problem of interest, we find that 
the induced current, obtained by Fourier transforming 
Eq. (Uni), is 

(r(q,c.))=™xf (q,^)u0- (14) 

xf is given by Eq. ^ with iq = and = 
where 

•^q = ^ E(2^" + 9a)4^Ck+q^ (15) 
ko- 

is the current operator with a denoting the different in- 
ternal states of interest (e.g., spin). 

Next, we need to relate Eq. to viscosity, using 
"constitutive relations" between the current and trans- 
port coefficients. For this we use Eqs. © and ([7]) sub- 
stituted into Eq. ([5]), where the symbols ja and Hap, 
without the hats used for operators, are understood to 
denote expectation values. In the long- wavelength limit, 
the contributions to the stress tensor coming from vis- 
cous terms dominate over contributions from pressure 
fluctuations, while the convective term dpUaUfj is be- 
yond linear order in velocity. We thus get mdtj" = 
(daCV ■u) + ri [V^-Uq -I- 9q,(V • m)/3] . Fourier transform- 
ing and comparing with Eq. (I14p . we obtain 

(16) 

We decompose the current correlation function into its 
longitudinal (xl) and transverse (xt) components: 

= + XT [5ap - J (17) 

By taking appropriate g — limits [24] of Eq. (fT6|) we 
find 

77(w) = lim {-iuj)m^XT{q, (18) 

q^O 

and 

C{uj) + 4?7(tj)/3 = lim {~iLo)m\L{q,io)/q^. (19) 

These expressions define the complex shear and bulk vis- 
cosities. We will be interested in the properties and sum 
rules of the spectral functions: 

Re ?7(w) = lim m^i:jImxT(q, w)/g^ (20) 



and 

ReC(w) + 4Re 77(a;)/3 lim m2wIniXL(q, a;)/^^. (21) 

g-i-O 

The static viscosities rj and C introduced in Eq. ([T]) are 
?7 = Re 77(0; = 0) and C = Re ({oj = 0). 

In closing this subsection, we note that the Kubo for- 
mulas for the viscosity derived here and below are valid 
in both the normal and superfluid phases, provided we 
recognize that the bulk viscosity in the superfluid state 
refers to ^2, which describes damping associated with 
an in-phase motion of the superfluid and normal fluid 
components [iTj . To understand this in more detail, we 
recall Landau's two-fiuid hydrodynamics JJJ for the su- 
perfluid state. In this theory, three bulk viscosities, (i, 
C2, and C.3, are required to describe the dissipation as- 
sociated with different types of relative motions of the 
superfiuid and normal components. The longitudinal re- 
sponse does not distinguish between the superfluid and 
normal components j22| and thus forces the superfluid 
and normal fluid velocities to be equal: Vg = v„ = u. 
When both components flow with the same velocity, the 
two-fluid hydrodynamic stress tensor [see Eq. (140.5) in 
Ref. [r?]] reduces to the expression in Eq. Q, with ^ 
replaced by C2, the bulk viscosity associated with the 
damping of the in-phase motions of the superfluid and 
normal fluid components. One can also show by direct 
application of Eq. (|82|) to the two-fluid hydrodynamic 
density response function in Eq. (|B2p that the left-hand 
side of Eq. is C2 +4/7/3 in the low-frequency two-fiuid 
hydrodynamic regime. 



B. Stress correlators 

We next derive Kubo formulas equivalent to those de- 
rived above but expressed in terms of the correlators of 
a suitably defined stress tensor operator Ha/s- These are 
useful to make connections with the literature [3, [2^ [2^ . 
We will also use these results in connection with the pos- 
itivity of the bulk viscosity spectral function and its van- 
ishing for the unitary Fermi gas. 

The IIq,^ operator must satisfy 

tm[ja,H]^dpflai3, (22) 

which is the operator version of the Euler equation, 
Eq. ([5]), and is simply a statement of momentum con- 
servation. We go to Fourier space and relate matrix ele- 
ments of the current operator to those of the stress tensor 
by sandwiching Eq. (|22p between exact many-body eigen- 
states. Using the spectral representation in Eq. we 
can then relate the current correlator x'j^i'i,'^) to the 
stress correlator Xn'''''''(qi '^)- The latter is defined by 
choosing A = n"'3(q) and B = m"^(-q) in Eq. (fT2)) . 

For simplicity we calculate only Xj^: which will suf- 
fice for our purposes. The final result, after some simple 
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algebra, is 



(23) 



Note that n' 



Q/3 



IIq^ + Aq^, with any symmetric ten- 
sor A satisfying dphap — 0, wih also be a solution to 
the Euler equation, Eq. ([2^ . This non- uniqueness in 
the definition of 11 does not affect our final results for 
the viscosity, related to Xj^i since a symmetric A with 
9/3 Aq^ = makes no contribution to Eq. (|23p . 

Using the decomposition given by Eq. (|17p . and taking 
the appropriate limits, we find 



rn^uj^ lim 



XT 



lim 



xy.xy 

Xn 



where we have taken and qz to zero before 



m Lo lim 



lim 



Xn 



-([nr,i%]) 



, (24) 
and 
, (25) 



where we have taken qy and q^ to zero before qx ■ We note 
that the commutators on the right hand sides of Eqs. ([M)) 
and ([25|) only affect the real parts of xt and x l and not 
the spectral functions of interest, shown in the next two 
equations. 

Using the Kubo formulas given by Eqs. ([20|) and (pT|) 
that were derived above, we find 



Eqs. ((20|) and ((2T|) . The two sets of equations are sim- 
ply related by the exact conservation law, Eq. (|22p . 
Above, we focused on the dissipative parts of the re- 
sponse, i.e., the real parts of the viscosities. Compar- 
ing Eqs. dill) and (HH) with Eqs. ^ and (US]), we see 
that the imaginary part of rj and the imaginary part of 
{Arj/3 + C) are not given by limi^_>o lini^^o Rexn^'^^/w 
and lim^_j.o limq^o RexJ^^'^^/w, respectively. Im r/ and 
Im C, when written in terms of stress correlators, also 
involve the frequency- independent, equal-time commu- 
tator terms in Eqs. and (|25p . This point seems to 
be missed in treatments that start out with the stress 
correlator formalism. The imaginary parts of the trans- 
port coefficients are most simply expressed in terms of 
the current correlation functions, Eqs. ^TE\\ and (IT^ . In 
the w — limit, the validity of this assertion can be 
seen quite independently from hydrodynamics (see Ap- 
pendix |B]). Allowing 1] to be complex in the hydrody- 
namic expression for the transverse current correlation 
function in Eq. (|B7p , for instance, one can readily confirm 
that the imaginary part of the shear viscosity is indeed 
given by Eq. p^ . 



III. SHEAR VISCOSITY AND NORMAL FLUID 
DENSITY 



Re rj{u}) = lini Imxn i^:^)/^ 



(26) 



and 



Re C(w) +4Re ?7(w)/3 = Ihn Imxn"'""(q,^)/^- (27) 

In an isotropic system, in the q limit, the only 
fourth rank tensor allowed by symmetry is of the form 
ASapSfj_i, + B {Safj.5pi, + Sav^Pti)- Wc cau thus combine 
Eqs. (IMl) and ^ to write 



Re C — ^Re r] 



Imxn'^^^''(q,w) 



lim 



(28) 



A very useful formula for the bulk viscosity follows 
from Eq. (|28p by looking at its {xx, yy) component and 
combining it with the {xx,xx) component in Eq. (1271) . 
Using the summation convention, we thus obtain 



Imxn"'^^(q,a;) 
9uj 



Re C('^) — lim 

g— >0 



(29) 



We emphasize again that the Kubo formulas for 
the bulk and shear viscosities expressed in terms of 
the stress-stress correlation function are equivalent to 
those expressed in terms of current-current correlations. 



In this Section we discuss the relation between the 
static shear viscosity Re 77(0; = 0) and the normal fluid 
density p„ , both of which can be written in terms of the 
transverse current-current correlation function. This al- 
lows us to prove that a non-zero normal fluid density p„ is 
a necessary condition for a non- vanishing shear viscosity 
77. This is, perhaps, not entirely unexpected on physical 
grounds, but we are unaware of a microscopic proof, valid 
for all Galilean invariant Bose or Fermi quantum fluids, 
that does not rely on a quasiparticle approximation. 

Before turning to the calculation, it may be useful to 
review some elementary facts about the shear viscosity 
77. Given that there is a Kubo formula for 77(0;) in terms 
of the current-current correlation function, Eq. (|20p . and 
that in kinetic theory rj is proportional to the mean free 
path, it may seem natural to assume that the shear vis- 
cosity of a fluid is the analog of metallic conductivity. 
This, however, is completely misleading. The shear vis- 
cosity is, in fact, the analog of the resistivity. This is 
clear, e.g., from the classical formula of Poiseuille for the 
flow rate Q = TrR^AP/iSrjL), with a pressure difference 
AP across a cylindrical pipe of radius R and length L. 
We will see below that zero viscosity in a superfluid is 
the analog of zero resistance in a superconductor. 

We begin by rewriting the Kubo formula for the shear 
viscosity, given by Eq. (j20p . using the spectral represen- 
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tation in Eq. (IT3|) : 



Keri(uj) — lim 

^ ' (T) Z 



Eh. 



a.b 



Siw-Eba), (30) 



Here and below, the "transverse limit", denoted by 
limj-T''), means that for Xj^ &rst set — and then 
take the limit qy ^ 0. 

The normal fluid density p„ characterizes the response 
of a fluid to moving walls and determines the moment 
of inertia of a cylinder containing the fluid; see, e.g., the 
detailed discussion in Refs. [12;[23- It is defined in terms 
of the real part of the static transverse current correlator: 



Pn — lim m Kexril,^ ~ 0). 



(31) 



Using the spectral representation in Eq. ([T2|) for Xj^, we 
can rewrite this result as 



-mfM\'- (32) 



^2 _ 

p„ = hm — - > — — 

(T) Z ^ Eta 



Our goal now is to understand the connection between 
the shear viscosity r/, which is obtained by taking the 
lim^_>o lim^_^o of Ito-Xt in Eq. (1301) . and the normal 
fluid density p„, which is the limq_>.o limt^^o of Rex^ 



in Eq. (|30|) . In lattice models of superconductors, it 
has been suggested [1^ that the order of the q and cj 
limits can be safely interchanged for the transverse cur- 
rent correlator, because all "transverse" excitations are 
gapped (unlike longitudinal excitations such as phonons 
in charge-neutral systems). However, this argument is 
not valid for the systems of interest to us. This can be 
seen, e.g., from the hydrodynamic form of xt in Eq. (jB7P 
which has a "diffusion pole" that makes the order of lim- 
its quite different. 

To prove the result stated at the beginning of this Sec- 
tion, we will show that pn = implies 77 = 0. The start- 
ing condition p„ — makes sense only at T = 0, since 
at any finite temperature there will necessarily be some 
thermal excitations. Furthermore, the vanishing of the 
normal fluid density 



Pr, — lim 2m 

^ (T) 



b 



EbQ 



(33) 



at T implies that each term in the sum over 
states vanishes. This means that, for each state \b), if 
the excitation energy varies as lim^y) _Ebo ~ with 
ab > 0, then the matrix element of the current operator 
vanishes even faster: lim(x) l(^bq|0)| ~ q°'b+Pb -v^jth > 
0. Note that we are not making any assumptions about 
the nature of the spectrum since both gapless (of, > 0) 
and gapped {at — 0) excitations are permitted. In either 
case, the matrix element of vanishes, since the q — > 



limit of jq is the total momentum, which commutes with 
the Hamiltonian in a Galilean invariant system. It is only 
in such a system that pn vanishes at T = [13, Ull ■ 

Now that we have constrained the matrix elements for 
any form of the excitation spectrum given pn — 0, we 
now ask how these constraints impact the shear viscosity. 
We look separately at the contribution from gapless and 
gapped states to Eq. (PU]) . which at T = can be written 
as 



-5{u-Ebo). (34) 



77(0;) — limTrm'^ Em- 



Each gapless state 6, with > 0, will contribute a term 
lim^ji-) (7^"''"'"^''~^i5(a; — A^q'^''), which gives a vanishing 
contribution [33| in the limit g — )■ for all w > 0. Finally 
taking the a; — > limit, we find that the contribution of 
the gapless states to 77 vanishes. 

Next, consider the gapped states with ab = 0, so 
that lim(3^) E^^ = A{, > 0. Their contribution to 
Eq. (p4| yields an expression of the form 77(0;) — 
lim(r) Y^'^^Cbq^''~'^5{ijj — At), where the prime indicates 
a sum over all gapped states. This result contributes to 
both the T] sum rule and the high-frequency tail that we 
will derive later in the paper. The important point here 
is that for < w < minj, {Ah}, i.e., below the minimum 
gap of all excitations, r]{uj) = 0. 

Thus, we conclude that the vanishing of the normal 
fluid density implies that the static limit of the shear 
viscosity vanishes as well: 77 = 0. This means that the 
Galilean invariant ground state of a superfluid has zero 
shear viscosity [31]. This is similar to the zero d.c. re- 
sistivity of a charged superconductor, as already men- 
tioned at beginning of this Section. There is, however, 
an important difference in that the vanishing resistivity 
persists all the way up to the transition temperature Tc- 
Even though there are normal fluid excitations in a su- 
perconductor, the infinite conductivity of the condensate 
"shorts out" the normal fluid in a superconductor. In 
marked contrast, in a neutral superfluid, 77 vanishes only 
at T = 0. For < T < Tc, even though a condensate 
exists, the normal fluid excitations give rise to a non-zero 
shear viscosity. 



IV. POSITIVITY OF SPECTRAL FUNCTIONS 



We simplify notation and write from now on 
77(0;) = Re ri{uj) and ({uj) = Re C(w), 



(35) 



unless explicitly stated otherwise. This should cause no 
confusion since we will not be dealing with the corre- 
sponding imaginary parts. Before deriving sum rules for 
77(0;) and C('^) in Section |Vl it is important to discuss 
here their positivity properties. Every time we say 'pos- 
itive' we actually mean 'non- negative', a term we find 
awkward for repeated use. We will show that 



77(w) > and ((uj) > VtJ. 



(36) 
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The simplest approach is to make expUcit use of the spec- 
tral representation. We will see that this is sufficient to 
prove the positivity of 77(0;), but not that of ({uj). To 
prove the latter, we will calculate the power absorbed 
by the fluid from an external velocity perturbation with 

Let us begin with Eqs. ([20)1 and (|2ip and use the spec- 
tral representation given by Eq. with Aq = and 

BJ^ — j^q. The transverse and longitudinal components 
are obtained, as usual, by taking suitable g lim- 
its HI. Using |(n|j^q|™)P > anda;[l-exp(-/3a;)] > 
for all w, we see that both culmxTil, w) and wlmxi(q, uj) 
are positive. Thus we obtain 



r]{io) > and C(w) + 4?7(w)/3 > Vw. 



(37) 



The inequality for C,{uj) is much weaker than what we 
wish to prove. One reason to expect that a stronger result 
should exist for C{uj) is that it is known from hydrody- 
namics (see Sec. 49 of Landau and Lifshitz [l3|) that the 
static bulk viscosity C(0) must be positive. To generalize 
this to all frequencies, we exploit the idea that the time- 
averaged power absorbed by the system from an external 
perturbation is necessarily positive. 

The rate at which the external velocity perturbation 
given by Eq. ([9]) does work on the fluid is given by 



dW 



dre-'"*e°"Vr)-(j(r,t)). (38) 



Following Ref . [32| , one finds that the time average of the 
power absorbed by the fluid is 



dW 
~dt 



Tfl 



q 



Q/3/ 



dW/dt > follows from the fact that energy can only 
be dissipated for any choice of the external velocity 
fleld. This implies that the real, symmetric matrix 
a;Im%"^(q, oj) must be positive definite, which is equiva- 
lent to the positivity of its eigenvalues. Using Eq. ([T7| . 
we see that these eigenvalues are precisely a;ImxL(q, t^) 
and wImxT(q, w), so that we simply rederive Eq. (1371) . 
and do not obtain ({uj) > 0. 

To constrain without any r]{uj) contribution, we 

must look at an external velocity field u{r, t) = M(r)e~*'^* 
with u{r) — ar, where a = (V • it) /3 is spatially uniform. 
To analyze the effect of such a perturbation, we first need 
to rewrite Eq. (|39|) in terms of the stress correlator so 
that daU/i is directly involved. Second, u{r) = ar is not 
Fourier transformable, so we must work in r-space, rather 
than q-space used elsewhere in the paper. 

We use the same derivation that led from the operator 
Euler equation given by Eq. ([^ to Eq. (^5]) . to get 



m^uj'^lmx"/{(l,uj) = q^g^lmxn^'''''(q,t^) 



(40) 



Using this in Eq. p9l) and rewriting the resulting expres- 
sion in real space, we get 



dW 
~di 



dr / dr'x 



Imxn^"^''(r - r'. 



dpu^ir'), (41) 



which must hold for arbitrary velocity fields u{r). 

To isolate the contribution of the bulk viscosity, we 
choose the velocity field u = ar, for which the shear 
term (in square brackets) in the viscous stress tensor, 
Eq. ([7]), vanishes. Using daUfj = aSap in Eq. (gT)) we 
get Imx^'^^ {<l — >■ 0,uj)/uj > 0. From the result given 
by Eq. ((29l) for the bulk viscosity, it immediately follows 
that C{uj) > for all cj. 



V. SUM RULES 

We now derive sum rules for the shear and bulk vis- 
cosities, dujri{uj) and dujC,{uj). We will first show 
that 



— / lim-^Imx^^(q 

([j-q,[^,iq]]: 



lim 

q->-Q 



2q2 



Then we will simplify the two terms on the right hand 
side of Eq. (|42|) : the first term by explicit evaluation of 
the commutators, and the second by appealing to hydro- 
dynamics. 

To see what is involved in deriving Eq. (|42|) . let us 
first be naive and ignore the q — >■ limit. Evaluating the 
integral on the left hand side by using the spectral repre- 
sentation in Eq. for Imxj^, we only obtain the first 
commutator term on the right. But taking the q — ^ 
limit after doing the cj-integration is not the same as in- 
terchanging the order of these operations! In order to 
do it correctly (g — > limit before the w-integration) , 
we exploit the Kramers-Kronig (K-K) relations to evalu- 
ate the integral in Eq. The only subtle point in this 
approach is that we need to ensure that the analytic func- 
tions which we K-K transform decay sufficiently rapidly 
for Lo — > 00. 

Using the expression in Eq. (|12|) . it is straightforward 
to expand the current correlator in powers of for 
large frequencies. One finds [33j . 



lim x7(q,( 



:[i-q,Jq]) ([i-qJ^,Jq]]> 



(43) 



The 



oj ^ term vanishes since ([jfiq,iq]) = 
~{'2qx/'m^)^]^cr n\^akx = in a uniform system. 
We further note that this expansion is strictly valid 
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only for a smooth potential [34| . a point which we will 
elaborate on in later Sections. 
Let us define a function F(uj) as 



F{uj) = lini — 

9^0 



xT 



, (44) 



where the q — > limit is defined appropriately [24| for 
the longitudinal and transverse cases. From Eq. (|43l) . we 
see that lim^^oo vanishes at least as fast as ui~^ 

and we can K-K transform it. We thus obtain 



lim ReF{uj) 



doj' 



,lmF{u') 



(45) 



dcj' lim ^ImY^T^(q, w') 



where we have used the fact that ajlmx^^(q, i^) is an even 
function of w. Using Eq. (|44p on the left-hand side of 
this expression immediately leads to the result, Eq. (|42l) . 
quoted above. 

As mentioned earlier, limi^^o linig^o('^^/2g^)Rexj^ in 
Eq. (j42p arises from the noncommutativity of the w — 
and g — >■ limits. Since this term involves the zero- 
frequency, long-wavelength limit where hydrodynamics 
is applicable, we can use hydrodynamic expressions for 
the current correlation function to evaluate it. In Ap- 
pendix [Bl we review such expressions and show that for 
any simple hydrodynamic liquid, one has 



lim lim 



2 2 

m ui 



RexT{q,uj) = 0, 



(46) 



lim lim 



2 2 

m uj 



9-)-0 



RexL(q,w) = 



where the adiabatic sound speed is Cg = (dP/dpY^^ at 
fixed s = S/N. Equation (|46)) is valid for both normal 
fluids and superfluids (within two- fluid hydrodynamics). 

Combining Eqs. dH]), and we flnd the 

following sum rules: 



1 



— duj 



lim 



^^([j-q,[g,iq]])T 

2g2 



(47) 



C(^) 



4ry(tj) 



: lim 



^^([AJg.J^]])L pel 

2q2 2 ■ 



Here, 



(48) 

)t(l) denotes the g — > limit appropriate to 
the transverse (longitudinal) case (23 |. 

The last remaining step in our derivation is to evaluate 
the commutators in Eqs. (j47|) and (j48|) . We consider a 
system of fermions or bosons described by the Hamilto- 



H ^ K + V 



(49) 



k+p<T'^k'-pc7'"^k'c 



■Cko 



k(T 



kk'p 



For a single-component Bose system, a = a' assumes one 
value; for fermions, a =t, i can take one of two "spin" 
values. It is straightforward, but tedious, to evaluate the 
commutator in Eq. for this Hamiltonian. One flnds, 
for both fermions and bosons. 



(50) 



8m 2 



- nip - q|)(p.-fe)' - nip + q|)(p.+&)'»- 

Here, {K) — X^ko- ^k'T^ko- is the kinetic energy density, 
and we have introduced the shorthand notation 



((2» - \ E 

kk'p 



.t 



k+po" k' — po" 



/Ck'fT'Ckcr). (51) 



Related expressions specific to Bose liquids are given in 
Ref. [111. We also note in passing that the longitudi- 
nal component of Eq. ([50)) is related by Eq. (jBll) to the 
so-called "(cj'^) sum rule" discussed for electronic sys- 
tems (36j . 

The right-hand side of Eq. (ISU]) varies as as g — 0, 
which cancels the 1/q^ in Eqs. (|47| and p8)) . Evaluating 
the transverse and longitudinal limits of Eq. ([50)) . one 
finds the following viscosity sum rules: 



1 



dijjrj(uj) 



3 



2V 
15" 



V 
30 



(52) 



and 



o 5e 
dwC(a;) = Y 



11^ + ^ + -^-^ (53) 
9+9+18 2 • ^^^^ 



Here, e = (ii') + (F) is the total energy density, (V^) the 

potential energy density, and the terms V and V are 
defined using Eq. (fSTj) as 



V = {{p{dV/dp))) and V = ((/ (9V/V))) • (54) 

These sum rules are valid at all temperatures (i.e., in 
the superfluid as well as normal phase) for any Bose 
or Fermi system with an arbitrary, spin-independent, 
isotropic interaction potential V{p). We emphasize 
that these are exact results obtained without making 
any quasiparticle approximations. In the next Section 
(Sec. IVI[) . we simplify these sum rules for the case of a 
two-component Fermi gas with short range interactions, 
which is of relevance to experiments on ultracold atomic 
Fermi gases with Feshbach scattering resonances. 

Before closing this Section, let us briefly discuss viscos- 
ity sum rules using the stress correlator representation. 
For the shear viscosity spectral function, the sum rule 



duirjiuj) 



= lim lim ^Rexri'^''''(q,w) (55) 
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follows trivially from the Kramers-Kronig relation. To 
show that this is the same as Eq. (|T7)) . we use Eq. (IMl) in 
the second line of Eq. (|55|) . One can rewrite the commu- 
tator in Eq. (|24l) using the Fourier transform of Eq. , 
and set w^Rexr/q^ to zero using the hydrodynamic re- 
suh, Eq. (gni), to obtain Eq. gT]). 



VI. DILUTE TWO-COMPONENT FERMI GAS 

We now specialize to the case of a two-component 
Fermi gas in the dilute limit, where the effective range tq 
of the potential (van der Waals at "long" distances, with 
ro ~ lOOflo) is much smaller than the s-wave scattering 
length a and the mean interparticle spacing kp^. (In 
typical experiments, kp^ ~ lfj,m and 500ao < \a\ < oo.) 
In the zero range limit rg — > 0, all physical observables 
are universal (ro-independent) functions of the energy 
scale ep (or length scale kp^) and the dimensionless pa- 
rameters T/ep (temperature) and l/ikpa) (interaction). 
We will show that for Fermi gases, the results given by 
Eqs. ([5^ and ([5^ of the previous Section, reduce to the 
simple expressions given by Eqs. ([2]) and ([3|) in the Intro- 
duction. 

Our main task is to calculate the terms V and V , 
involving gradients of the interaction potential, defined 
in Eq. (|54)) . We use the real-space approach developed by 
Zhang and Leggett [37], which is a simple way to derive 
results first obtained by Tan [13, [11]. Using the two-body 
density matrix 

J-(r) = 



d^R(V|(R+^)V'l(R-x)V'4(R-;T)^t(R+;j) 



we rewrite V and V in real space as 

-J-(r) 



T7' I j3 ^^(0 

V = (TV r- 



and 



V = d'v 



dr 



Qj.2 



(56) 



(57) 



(58) 



Since V{r) is short-ranged, these expressions are only 
sensitive to the short-distance (rp < r <^ kp^) structure 
of the two-body density matrix. (The non-universal con- 
tribution from distances smaller than is assumed to be 
small.) For a two-component dilute Fermi gas, at these 
short distances, the two-body density matrix is (37| 



.F(r) 



C 
167r2 



(59) 



Here, C is the contact [T3, 113, HI] mentioned in the In- 
troduction. In Appendix [C] we remind the reader how 



the contact C governs both the short-distance behav- 
ior of the two-body density matrix in Eq. (|59p . and 
the large-fc tail of the momentum distribution function 
limfe^oo "kcr ^ C/k^. 

Using integration by parts, we transform gradients of 
the potential V{r) in Eqs. (|57|) and ([58]) into gradients of 
the two-body density matrix, Eq. ([59| . We thus find 



(60) 



and 



V' = — f drV(r) (-1 + 4r/a) 
47r J 



V" I drV{r) (1 - 6r/a) 

ZTT 



(61) 



All that remains is to evaluate the two integrals X„ = 
C J drV{r){r/a)'^ /Att with n = 0, 1 in the limit where the 
range of the potential rg — ?> 0. The Tan relations are pre- 
cisely what we need to evaluate such (possibly divergent) 
integrals. The details of this analysis are described in Ap- 
pendix [Cl We use the potential energy density [13, HI] 



(V) 



c 



CK 



(62) 



where A = 1 /rg is the ultraviolet cutoff , and the pressure 

P = 2£/3 + C/(127rma) (63) 

to determine Xq and Xi. In deriving these results, we 
also use an expression for the pressure P in terms of e, 
(V) and V which is derived in Appendix |D] using the 
Feynman-Hellmann theorem. 

Our final results for V and V , derived in Appendix 
O are 



V = -{V) -2e + 3P = Ckj^Tp-m 



and 



V = 2(1/) +8e- 12P = - 



CA 



C 



2TTma 



(64) 



(65) 



Using these results in the general sum rules given by 
Eqs. (|47p and we obtain the 77 and C sum rules for 
the two-component dilute Fermi gas which are valid for 
all values oil/^kpa) throughout the BCS-BEC crossover, 
so long as a, kp^ 3> rg, and at all temperatures, both in 
the superfluid and normal phases, so long as T <C 1/rnrg. 
For the shear viscosity, we find 



dujri{uj)/Ti = e/3-2(y)/5 



C 



CA 



(66) 



where we have imposed the energy cutoff A^/m = 
1/mrg [1^. In the zero-range limit as A = l/rg ^ 00, 
the right hand side diverges. (Strictly speaking, every 
physical potential has a small non-zero effective range rg. 
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which leads to a weh-defined, finite results, but one that 
is "non-universal" in that it depends on short distance 
physics.) We will see in the following Section, Sec. IVIH 
how to make sense of this divergence and find a modified 
sum rule that remains finite as tq — >■ 0. 
For the bulk viscosity we find 



dwC(w)/7r = P - e/Q- pcl/2 



5£ 



c pel 
Ylnma 2 



(67) 
(68) 



Below the superfluid transition, the bulk viscosity C that 
enters Eq. ([3]) is the bulk viscosity (^2, as explained earlier. 
We can rewrite the right hand side of this sum rule in 
a useful way using simple facts about the scaling form 
of thermodynamic functions across the entire BCS-BEC 
crossover, as described in detail in Appendix|El The final 
result is 



duj C,{Lo)/Tr 



1 



/ dC 



(69) 



where the derivative is taken at constant entropy per par- 
ticle s = S/N . The positivity of the sum rule, given that 
of its integrand, implies that the contact is a monotoni- 
cally increasing function of l/{kpa) through the BCS- 
BEC crossover. We will discuss this further in Sec- 
tion Eml 



VII. HIGH-FREQUENCY TAILS 

In this Section we derive a modified shear viscosity 
sum rule that is manifestly finite in the A — l/rg — >■ 
oo limit. This is obtained by relating the linear (in A) 
divergence in the sum rule, Eq. (|66p. to a high-frequency 
tail in rj{ijj) ^ 1/v^i and then "subtracting out" the 
contribution of this tail. We use "high frequency" or 
w — >■ oo, to mean ei? ^ a; < 1/rnrg. We also argue that 
a high-frequency tail of the form w^"/^, with odd integer 
n, in a variety of spectral functions is a generic feature 
of short-range physics. As discussed below, it shows up 
in many contexts, even outside dilute quantum gases. 

We can rewrite the rj sum rule in Eq. (j66|) as 




lOTTJTia 



(70) 



where fig is an arbitrary energy scale. If we choose i7o 
to be A? /m we recover Eq. (|66p . But for any finite 
subtracting out the lo~^/'^ tail makes the integral ultravi- 
olet convergent and we can take the cutoff A to infinity. 
If we choose fJo = 0, we obtain the finite sum rule 



— I duj 



C 



C 



IOtttoo 



(71) 



The price we pay for using this finite, ro-independent re- 
sult (in the rg — >■ limit) is that we sacrifice the positivity 
of the integrand. At sufficiently small w, we must nec- 
essarily have ri{Lu) < C/ {lOTTy/niLu) since 77(0) is finite. 
One can, in principle, exploit the freedom in Eq. (j70p 
and choose flo to be large enough so that the integrand 
is always positive, however. 

The finiteness of the right hand side of Eq. ([7T|) implies 
that the integrand on the left must vanish at least as fast 
as for the integral to converge at large uj. Thus 

the asymptotic behavior of the spectral function 77(0;) is 
of the form 



77(0-1 — > 00) 



C 



107^^ 



(72) 



We note that a high-frequency tail in the imaginary 
part of a retarded correlation function which goes like 
with positive integer n, is a general feature of 
short-range two-body physics. Suppose that for some 
operator A, the corresponding 7i-th moment sum rule 
has the form 



1 



dLULL>^lmxA,Ai^) — ce{V) + 



(73) 



where we only show the divergent term explicitly; the 
ellipses denote regular terms, a is some combination of 
parameters and is not, in general, dimensionless. In addi- 
tion to the current correlation function (n = 1), diverging 
sum rules of the form given by Eq. ([73]) arise for the ra- 
dio frequency (RE) spectral function {n = 1) [40], and, 
as we show below, the density response function (n = 3). 
Using the same reasoning as above, a divergence of the 
form given by Eq. (I73p implies a high-frequency tail. For 
a dilute two-component Fermi gas with a 3> t-q , the high- 
frequency tail is given by 



ImxA.Aii^ -> 00) 



1 



(74) 



As seen from the above arguments, an uj~^^^ tail arises 
in the radio-frequency spectroscopy response function 
I{uj) [41I, \^ for Fermi gases. Another interesting ex- 
ample is the w"^/^ tail in the density response of a dilute 
Fermi gas which we derive in Section [IXI There, we also 
point out that an identical asymptotic behavior is found 
for the dense Bose liquid ''He, which further emphasizes 
the generality of the short-distance physics in all quan- 
tum fluids. 



VIII. SUM RULES THROUGH THE BCS-BEC 
CROSSOVER 

In this Section we consider the bulk and shear vis- 
cosity sum rules through the BCS-BEC crossover, going 
from the weakly attractive BCS limit (a small and nega- 
tive) with large Cooper pairs to the BEC limit (a small 
and positive) with weakly interacting, tightly bound 
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FIG. 1: (Color online) The value of the bulk viscosity sum 
rule given by the right-hand side of Eq. (1691) at T = in units 
of nep through the BCS-BEC crossover. 



molecules. The crossover can be traversed by changing 
X = l/{kpa) from x — — oo (BCS limit) to x = +cxi 
(BEC limit). In experiments, the scattering length a is 
varied by tuning a magnetic field about a Feshbach res- 
onance. Precisely at resonance, a; = 0, the scattering 
length diverges and the Fermi gas is in a very strongly 
interacting "unitary regime" where the pair size is of the 
order of the interparticle spacing. 

To actually compute the viscosity sum rules given by 
Eqs. ([M]) and ([7T]) for arbitrary coupling x — l/{kpa) 
and temperature T, we need to know the energy den- 
sity e — neF£{x,T / ep), from which we can determine 
the contact C as described below. In general, we will 
need to use quantum Monte Carlo (QMC) data for the 
energy density to evaluate the sum rules. However, as 
shown below, we are able to analytically constrain the 
bulk viscosity spectral function at unitarity. 

We see from Eq. that the bulk viscosity sum rule 
vanishes at unitarity: 

1 

- I duj C(w) ^ (|a| = oo). (75) 
^ Jo 

We are using here the fact that {dC/da~^)s is finite (i.e., 
non-infinite) at x = l/^kpci) = at all temperatures. 
One can also see this using elementary arguments that 
do not involve the contact. From "universal thermo- 
dynamics" [i^, the only energy scales at unitarity are 
ep and the temperature, and we can directly show that 
P-e/9- pcl/2 = (see Appendix |E|. 

The vanishing sum rule, Eq. (|75p . together with the 
positivity condition C(^) ^ derived in Section llVl im- 
plies 

C(w)=0 Vw {\a\=oo). (76) 

That the static bulk viscosity (^(0) vanishes is a well- 
known consequence [l5| of scale or conformal invariance 



0.4 




FIG. 2: (Color online) The value of the finite shear viscosity 
sum rule, with the contribution from the high-frequency tail 
in Eq. (|72|l subtracted out, given by the right-hand side of 
Eq. ^ at r = in units of nep through the BCS-BEC 
crossover. 



at unitarity j44i |. Our result generalizes this to arbitrary 
frequencies. As discussed below in Section HXl our result 
actually has important implications for measuring the 
frequency dependent shear viscosity of a unitary Fermi 
gas using a density probe such as two-photon Bragg scat- 
tering. 

Another general consequence of ({lu) > is that its 
sum rule must be positive for all x = l/{kpa) and T. 
Equation then implies that 

so that the contact must be a monotonically increasing 
function of l/{kpa) through the BCS-BEC crossover at 
fixed entropy per particle. We can understand this in- 
equality intuitively as follows: the contact C, which is 
related to the probability of finding two particles of op- 
posite spin close to each other, can only increase with 
increasing attraction a"^. 

In Fig. [1] we show the bulk viscosity sum rule in 
Eq. dnH) at T = calculated using QMC data 'IS] for 
the energy density e. The contact C is obtained from e 
using Tan's "adiabatic relation" [l^ 

{de/da-^)^ = -C/(47rm), (78) 

where the derivative is taken at fixed entropy per particle 
s = S/N. We fitted the QMC data and took numerical 
derivatives with respect to a~^. Since the ^ sum rule 
involves the second derivative of QMC data for the energy 
density, the results may not be very accurate far from 
unitarity in either direction. 

Both the vanishing of the C sum rule at a; = l/{kpa) — 
and its positivity away from unitarity are apparent in 
Fig. [1] This is due to the 1 /a^ dependence of the sum 
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rule in the vicinity of unitarity. We emphasize the non- 
triviahty of the result given by Eq. ([S^ in the unitarity 
region. In the form first derived in Eq. (|68|) . the right 
hand side is (P — e/9 — /9c^/2). Each term in this expres- 
sion has both constant and order x contributions, which 
must all cancel to give a final result which goes like a;^ 
at small x. In the BCS limit, the C sum rule vanishes as 
2neF/i27TT\x\) since C iir^n^al [3. In the BEC limit, 
the energy density is dominated by the negative molec- 
ular binding energy, e « nEi,/2, with E}, = —l/{mo?). 
Thus, C — T' Airn/a and the sum rule grows as n\Ei,\/l%. 

Next, in Fig. [51 we plot the shear viscosity sum rule 
given b y E g. (|7T|) at T = again using the QMC data 
of Ref. [43 ■ Because of the 1 / subtraction extending 
all the way down to w = 0, the 77 sum rule in Eq. (j7ip 
is not constrained to be positive. Using the above an- 
alytic result for the contact in the BCS limit, one finds 
that the 77 sum rule asymptotes to Q.2nep in the BCS 
limit. At unitarity, \a\ — 00 and the 77 sum rule is 
e/3 ~ 0.4 X {3neF/5) x (1/3) = 0.08^6^. On the BEC 
side of the resonance the sum rule changes sign, tending 
to (17/30)n£'f, in the BEC limit. 



IX. DYNAMIC STRUCTURE FACTOR 

We now discuss the connection between viscosity and 
the density-density correlator or dynamic structure fac- 
tor. This analysis leads to two interesting results for the 
two-component Fermi gas. First, we predict that a den- 
sity probe such as two-photon Bragg spectroscopy 
can in principle be used to measure the frequency depen- 
dent ri{uj) at unitarity: 

T]{uj) = \im — Imxpp(q,w) (|a| = 00). (79) 

Second, we derive the high-frequency tail [l6j 

2q'^C 1 



lim lim ^(q, oj) = — 

uj^oo f/— >0 15 



(80) 



a result that is valid for all l/{kFa) and all tempera- 
tures. As discussed below, such non-analytic tails are 
also known in other strongly interacting quantum fluids 
like ''He. 

We start with the operator form of the continuity equa- 
tion 



i[p,H] = mdaja 



where p — mh is the mass density operator, and take its 
matrix elements between exact many-body eigenstates. 
Using the spectral representation, Eq. (fT2l) . we relate 
the density correlator Xpp to the the longitudinal cur- 
rent correlator [see Eq. (jBl[) ]. The latter is related to the 
viscosity as shown in Eq. ([2T|). namely ({u) + 4:r]{uj) /3 — 
liniq^o ?'T'^'^IniXL(q, w)/g^. We thus obtain 

C(c.) + 47,(c.)/3 = lim io^lmXpM^ cu)/q\ (82) 



We discuss two situations where the contribution of ({uj) 
vanishes and we can obtain interesting results connecting 
ri{uj) and density correlations. 

First, we focus on the unitary Fermi gas where Ci'-^) 
vanishes at all u! (as shown in Section IVlIII) and Eq. ([5^ 
simplifies to Eq. (1791) . Thus, the frequency-dependent 
shear viscosity i]{uj) in a unitary Fermi gas can in princi- 
ple be measured using an experiment like Bragg scatter- 
ing, which directly probes Imxpp. 

Second, let us look at the high-frequency regime ej? <C 
a; < 1/ (ynrg). The C sum rule in Eq. (I68|) is convergent in 
the ro — >■ limit, and thus C{uj) must decay faster than 
1/uj, while ri{ijj) ~ l/Vw (see Section IVIl|) . Thus, as 
cj — > 00, the bulk viscosity Ci^) is much smaller than the 
shear viscosity ri{bj) for all l/{kpa) and all T/ep. Using 
Eq. dHU) we thus find 



T]{iU — > 00) 



3^3 

w— S-cxD q— >o Aq 



lini_ hni ^j-jlmxpp(q, uj) 
37rcj3 



lim lim 



-5(q,cj), 



(83) 



The dynamic structure factor 5(q, oj) is related to Im^pp 
via the fluctuation-dissipation theorem 



IinXpp(q,^) 
7r[l — exp(— /?w)] ' 



(84) 



Our final result, given by Eq. (j80p . for the high-frequency 
tail of S{q,uj) is obtained by using the high-frequency 
tail, Eq. of 77(0;) in Eq. 

The high-frequency w"^/^ tail of the dynamic structure 
factor result is a universal feature of short-range two- 
body interactions. Remarkably, such a tail was first no- 
ticed in deep inelastic neutron scattering studies of super- 
fluid ■'He [431 and was subsequently understood in terms 
of hard-sphere gases [i^. The high-frequency neutron 
scattering experiments probe the short distance proper- 
ties of the two-body pair distribution function. [In dilute 
Fermi gases, this is directly related to the contact C; see 
Eq. (j59p ]. It may seem surprising that such anomalous 
high-frequency tails arise even in dense systems like *He. 
Recall that this behavior should be visible in a frequency 
range n'^^^/m < uj < l/mrp in which the interaction 
"looks" short-range. Even in "'He, where nr^ < 1, such 
a frequency range can be found using deep inelastic neu- 
tron scattering, although the range is obviously much 
smaller than in dilute gases with uvq ^ 1. 



X. COMPARISON WITH SUM RULES FOR 
RELATIVISTIC FIELD THEORIES 

There has been a considerable effort in the high-energy 
literature to understand the properties of viscosity spec- 
tral functions and their sum rules; see, e.g., Refs. |49l - 
Isijl). In addition to understanding the transport coeffi- 
cients within the AdS/CFT framework, this work seems 
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to be motivated in part by an interest in reliably ex- 
tracting transport coefHcients of the quark-gluon plasma 
from lattice QCD calculations of Euclidean correlation 
functions. We briefly discuss here some similarities and 
differences between the results for relativistic quantum 
field theories and those derived in this paper for non- 
relativistic Fermi gases: Eqs. ([M)) and (|71l) . 

There exist a number of Boltzmann calculations of the 
viscosity spectral functions in weak couphng QCD j2l.l50|. 
For the shear viscosity, the authors of Ref. 2] find the 
shear viscosity sum rule 

1 f^" e + P 

- / dujiico) - (85) 
Jo 5 

where g'^T ^ Uc g^T is a cutoff that removes a di- 
verging contribution from a high-frequency tail. For the 
jV = 4 supersymmetric Yang- Mills theory (SYM), Ro- 
matschke and Son [i^ derived the following shear vis- 
cosity sum rule: 

1 f°° e 

- duj[ri{uj) - riT=o{i^)] = -r- (86) 
Jo 5 

Here, a diverging vacuum contribution from a T- 
independent high-frequency tail has been subtracted out. 
We note that our 77 sum rule in Eq. ([7T|) , though similar in 
structure, has one key difference. The high-frequency tail 
for the Fermi gas is in general T- dependent, because its 
coefficient is set by the contact C = kpC[l/{kpa),T/eF]- 
A non-perturbative calculation of the bulk viscosity 
sum rule in A/" = 4 supersymmteric Yang-Mills theory 
and pure Yang-Mills theory (QCD with no quarks) has 
been given recently by Romatschke and Son [49i] : 

1 f°° 

- / duj[C{uj) - Ct=o(w)] = 

7^ Jo 

(3e + P)(l-3c2)-4(e-3P), (87) 

where c = ^JdPjde is the sound speed in relativistic hy- 
drodynamics (with the speed of light equal to unity) [13] ■ 
There are some differences and one very interesting simi- 
larity with our C sum rule in Eq. (|69p. In contrast to the 
Fermi gas spectral function C^iuS)., there is a need to sub- 
tract out a divergent tail in Eq. |57] and this tail appears 
to be T-independent. The interesting similarity is that 
in the "conformal limit" P = e/3, the right hand side of 
Eq. ([57)1 vanishes, analogous to the unitary Fermi gas. 

XI. CONCLUSIONS 

In this paper we have derived various exact, non- 
perturbative results for the shear and bulk viscosities of 
non-relativistic quantum fluids, focusing on the strongly 
interacting Fermi gas. Our main results were already 
summarized in the Introduction. To conclude, we dis- 
cuss some open questions and how our results relate to 
them. 



Most calculations [2^ [S^ of the viscosity in strongly 
interacting Fermi gases have so far been restricted to solv- 
ing Boltzmann equations or using diagrammatic pertur- 
bation theory, in essence making a quasiparticle approx- 
imation. Such results are valid in the high and low tem- 
perature regimes, but not in the most interesting regime 
near and above where a quasiparticle approximation 
is questionable and the shear viscosity is known to be 
the smallest. It was recognized some time back [sl] that 
there is a breakdown of Fermi liquid theory in the nor- 
mal (i.e., non-superfluid) state of the strongly interact- 
ing regime of the BCS-BEC crossover. It was shown 
that precursor pairing correlations lead to a pseudogap 
[53I I , which is a strong suppression of low-energy spectral 
weight in various response functions. It is likely that no 
sharp quasiparticle excitations exist in this regime near 
unitarity and just above Tc, but controlled calculations 
of dynamic quantities are very difficult. 

Quantum Monte Carlo methods have played an impor- 
tant role in determining the equilibrium thermodynamic 
properties of the unitary Fermi gas. However, results 
for transport coefficients are much less common, since 
they require analytic continuation of imaginary time (Eu- 
clidean) data to the real axis |5j|. The sum rules we 
derive could serve as useful constraints on similar calcu- 
lations for strongly interacting Fermi gases. 

From an experimental point of view, the (static) shear 
viscosity for strongly interacting Fermi gases has been 
estimated from studies of the damping of collective oscil- 
lations |11|]. We have shown above that, at unitarity, the 
full frequency dependence of the shear viscosity spectral 
function r]{uj) can be obtained from two-photon Bragg 
spectroscopy. While it would be a challenging experiment 
(the density response being very small for small-g), this 
would give extremely important insights into the strongly 
interacting Fermi gas, analogous to optical conductivity 
measurements of solids. 

Finally, we return to the conjectured bound jj] on the 
shear viscosity. Eg. ([I]). Proving or disproving the exis- 
tence of a bound 55] for non-relativistic quantum fluids 
like the strongly interacting Fermi gas remains a challeng- 
ing open problem. We hope that the spectral functions 
and sum rules derived here constitute a step in this direc- 
tion, just as they have for other well known inequalities 
in quantum many-body physics. 
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Appendix A: Modified stress correlator 

In the main paper we discussed current correlator and 
stress correlator representations of the bulk and shear 
viscosities. Here, we describe a third correlation function 
using an explicitly defined operator IlJ^^ which has been 
used to calculate the static shear viscosity rj — Ke t]{uj — 
0). We note that H", which does not include the diagonal 
terms of the full stress tensor, cannot be used to calculate 
the bulk viscosity. 

Let us define 



n°.(r) 



N 



i=l 



(Al) 



where is the a-component of the momentum opera- 
tor for the i-th particle. We emphasize that this is only 
one piece - the kinetic part - of the full stress tensor 
operator, and omits other terms, such as the pressure. 
It is independent of the interaction potential unlike the 
full stress tensor. However, since the expectation value 
of the "off-diagonal" part of H° is identical to the hydro- 
dynamic stress tensor in Eq. ([B]), we expect that we can 
use H° to compute the shear viscosity, at least in the low 
frequency limit. 

We define the correlator Xuo'^^ by choosing A — 
n°^(q) and B = fl^i-q) in Eq. m, where 

^%i<l) = ^ E cLck+q.(2/c„ + q^)i2kp + qp). (A2) 

kcr 



We can write the analog of Eq. as 
?7o(w) = lini Imx^^o'^^^/w. 

g— >0 



(A3) 



This is the form used by several authors |25|, |26[ as a 
starting point for diagrammatic approximations. 

The sum rule for the modified stress correlator, 
Eq. (jA2l) . simply follows from the Kramers-Kronig re- 
lation: 



lim lim — 

bj^o (j->o 2 



Rexno'''^(q,w) 



(A4) 



Ironically, it is seems harder to explicitly evaluate the 
right hand side here than it is to calculate the exact sum 
rule in Eq. (1551) . despite the simple operator H° involved. 
The point is that fl*^ does not satisfy the Euler equation, 
Eq. (|22|) . and hence we cannot relate it to the current. 
Thus, in contrast to the sum rules given by Eqs. (j47p 
and (|48p which involve the first frequency moment of the 



current correlator, we must deal directly with an inverse 
frequency moment in Eq. (|A4p . Such an inverse moment 
is a generalized "static susceptibility", about which we 
do not seem to know much, at least in this case. Unlike 
positive moment sum rules, it cannot be written in terms 
of commutators. 



Appendix B: Hydrodynamics 

In this Appendix, we review well-known hydrodynamic 
results for the current correlation functions [l^, [s^ . To 
keep the discussion as general as possible, we will use the 
full two-fiuid hydrodynamic correlation functions that re- 
sult from solving the linearized equations of two-fiuid hy- 
drodynamics p7l.[57|. As written below, these correlation 
functions describe any superfluid with a two-component 
order parameter, including dilute two-component Fermi 
gases (see Ref. [53] and references therein), and reduce to 
standard hydrodynamic expressions in the normal phase 
above Tc- We start by writing down a relation between 
the longitudinal current correlation function Xi(q, oj) 
and the (mass) density response function Xpp(q, w). (Re- 
call that our current correlation function is the number 
current correlation function and not the mass current 
correlation function generally used in the older litera- 
ture [1^ [5^ . We will find it convenient in the analysis 
below to use the correlation function Xpp fo'' the mass 
density p — mn, however.) Analogous to the result given 
by Eq. (^5]) . the continuity equation dtn -|- V • j = can 
be used to find 

^'^Xppi<l,^) = (»7tq)^XL(q,w) - q- ([jq,p_q]) 

= {mq)^XL{(l,(^)- pq^- (Bl) 

We can now use the hydrodynamic expression for 
Xpp (q,a;) to obtain an explicit hydrodynamic expression 
ioT XL(q, w). In the hydrodynamic regime, the density 
response function is (see, e.g., Eq. (4.32) in Ref. [56j): 



Xpp(q,'^) 



pq 



{uj'^ — ulq'^ + iDiq'^uj){uj'^ — u^q^ + iD2q'^uj) ' 



(B2) 



Here, ui and U2 are the speeds of first and second sound, 
respectively. They can be shown to satisfy the following 
identity (see, e.g., Eq. (14.39) in Ref. 

ul + ul = ^^ + cl with cl = {dP/dp)^. (B3) 

Pn^v 

RecaU that s = S/N is the entropy per particle. Cv = 
T[ds/dT)p is the specific heat per unit mass at constant 
volume, ps and p„ are the superfluid and normal fluid 
densities. The damping coefficients F, and D2 obey 
the following identities: 



K 

pcv 



Ps fAr,/5 + C2 



Pn 



P 



- 2Ci + PCs 



(B4) 
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and 



Z?i+i?2 = r+-(4?7/3 + C2). 

p 



(B5) 



Here, k is the thermal conductivity and Ci, C2j and Cs are 
the bulk viscosities associated with the different types of 
motion that can arise in the superfluid phase [l3| • Above 
Tc, C2 reduces to the usual bulk viscosity C, and the re- 
maining bulk viscosities do not contribute. 

After some straightforward but lengthy algebra, one 
can show from Eqs. (jBl|) and (|B2|) that 



--- ----- o -J-^^^L (q, w) = - ^ 



lim lim 



Pn 



(B6) 



Using Eq. (|B3| in this expression gives the result in 
Eq. dUl). 

The transverse current correlation function is given by 
(see Eq. (4.49) in Ref. ^) 



m 



(B7) 



From Eq. (jB7l) , we see that the real part of the transverse 
current correlation function is proportional to g**, leading 
to the result in Eq. 



Appendix C: Contact 

This Appendix consists of two parts. In the first part, 
we briefly recall, for completeness, some basic properties 
of Tan's contact; more details may be found in the orig- 
inal references [3 HI]. In the second part, we use the 
same techniques, within a real-space formulation [37|, to 
derive Eqs. (|64|) and (|65)) for v' and v" . 

The contact C can be defined by the large-fc tail 
of the momentum distribution function limj;_j.oo f^k — 
C/k^. This leads to a kinetic energy density {K) = 
2 / d'^k (fc^/2m) with a linearly divergent piece that 
goes like CA/(27r^m), where A = l/rg is the ultraviolet 
cutoff. The potential energy density is given by 38] 



{V) = 



c 



AC 



(CI) 



so that the total energy density e — (K) + {V) is finite 
in the A = l/rg — > 00 limit. We will freely use these 
results, together with those obtained from the short- 
distance properties of two-body density matrix, to eval- 
uate the quantities of interest for our sum rules. 

At short distances, tq < r ^ kp^, the two-body den- 
sity matrixfor a two-component dilute Fermi gas has the 
structure 37] 



)ViI(R-|)Vi;(R-|)V3t(R- 



^ J-(r) = 



C 
16^ 



---) -(02) 
r a ' 



Let us use this to compute the interaction energy density 

= drV{r){l-r/af. (C3) 

It is easy to see that for rp — >■ 0, we may drop the [r/aY 
term in the integrand as it gives a vanishingly small con- 
tribution. Using 

^" " £ / '^''^^''^ ^^Z''^" = 0, 1) (C4) 
we thus obtain 

(T>)=Ao-2Ai. (C5) 
Similarly, Eqs. (j60|) and (|6ip may be written as 

V' = -Xq + AXi (C6) 



and 



V =2X0-12X1. 



(C7) 



Next, we wish to determine the integrals Aq and Xi 
in the limit where the range of the potential vanishes: 
ro — > 0. Comparing the results given by Eqs. (|C1|) and 
(IC5[) for {V), it is evident that Aq is linearly divergent in 
A and Xi is finite as A = l/rp — > 00. But there is no way 
to determine the finite part of Xq from this comparison 
alone. We need one additional piece of information to 
determine Xq and Xi. We get this from the relation 



P = 2e/i+{V)/i + V /3. 



(C8) 



for the pressure P which is derived in Appendix |D] Using 
Eqs. (jC5|) and (jC6l) . we see that the divergent terms (Aq) 
cancel and the pressure 



P = 2e/3 + 2Ai/3 



(C9) 



is finite as — ^ 0. Comparing this with [TJ] P = 2£/3 - 
C / {\2ti ma) we obtain 



Xi = C/(87rma). 



(CIO) 



Substituting this into Eq. (jC5|) for (V^), and comparing 
with Eq. we find 



Ao = -CA/(27r2m) - C/(27rma). (Cll) 



Using these results for Aq and Ai in Eqs. ()C6|) and (|C7P 
we obtain Eqs. (|64l) and (|65|) for V and V , respectively. 



Appendix D: Pressure 

In the main text of this paper, we have suppressed 
factors of the volume at all places by setting it 
equal to unity. In this Appendix, we re-introduce fac- 
tors of fl in order to use the thermodynamic relation 
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P = —{dF/dVt)T,N to derive a microscopic expression 
for the pressure. To evaluate this, we use the Feynman- 
Hellmann formula {dF/dX)T,N = {dH/dX) treating the 
volume as a parameter A in the Hamiltonian 



H 



1 



Ck<T + 7: > , Vip)c- 



kk'p 



k+pt'-k' 



,p4Ck'4Ckt(Dl) 



The volume enters in two ways: (i) explicitly, though 
the r2~^ factor in front of the interaction term, and 
(ii) implicitly, through the discrete wavevectors k — 
(27r/r2^/^) {nx,ny,nz). One can replace the momentum 
sums by sums over the discrete indices tIq, = 0, 1, and 
the operators only depend on these indices. Thus, in 
addition to the explicit factor, the kinetic energy 
Ek = k^/2m and the two-body potential V{p) also de- 
pend on the volume fi. 

Using 3n[dF{k)/dn] = -ko,[dF{k)/dka] (with the 
summation convention) and the definition of V given in 
Eq. (|54|) . we find the result 



2 1,-, 1— ' 



(D2) 



Appendix E: Thermodynamics of the BCS-BEC 
Crossover 



s = S/N. We find it convenient to use s, rather then 
more familiar variable T/ep, because we will need to 
evaluate adiabatic derivatives below. The density de- 
pendence of the Fermi energy and Fermi wavevector are 
given by ep — (37r^n)^/'^/2r7i and kp — (STr^n)^/"^ respec- 
tively, and a is the s-wave scattering length. 

We first calculate the adiabatic sound speed which en- 
ters the right hand side of the C sum rule in Eq. (|68l) . 
Using the definition cj. = {dP/dp)s, where p = ran, to- 
gether with Tan's pressure relation, Eq. (1551) . we find 



PCs 

2 



dn 



1 



dC\ 

127rma V dn ) 



(E2) 



The derivatives at constant s are evaluated as follows. 
The first term is {de/dn)s = {5e/3n) — {epx/3){d£ /dx)s- 
We compute {d£/dx)s using Tan's adiabatic relation, 
Eq. dZHl), and obtain n{de/dn)s = 5e/3 - C/{12nma). 
To calculate the second term in Eq. (|E2[) . we rewrite the 
contact in the scaling form 



C = k%Cix,s) 



(E3) 



where C is a dimensionless function on its arguments. 
After some simple algebra, we find n{dC/dn)s = 4C/3 — 
{3a)^^{dC /da~^)s- Adding up all of the contributions 
to the bulk viscosity sum rule, we find 



In this Appendix, we simplify the form of the bulk vis- 
cosity sum rule in Eq. (|68l) and derive the result given 
by Eq. (|69p using thermodynamic scaling arguments [i^ 
and the Tan relations [l3|- We begin by writing the en- 
ergy density of a two-component Fermi gas in the scaling 
form 



£ = nep£{x, s). 



(El) 



which is valid across the entire BCS-BEC crossover. 
Here, f is a dimensionless function of the interaction 
parameter x = l/{kpa) and the entropy per particle 



he 



c 



12iima 



2 



1 



dC 



72'Kma?' \da ^ 



(E4) 



That the C sum rule vanishes at unitarity can also 
be seen directly from thermodynamic scaling arguments, 
without introducing the contact. Using Eq. (jEip . we see 
that P — —{d{en)/dV,)s N = 2e/3 at unitarity [as an- 
ticipated by Eq. (|63)) ] |59| . Using this, we also find that 
the adiabatic sound speed at all temperatures is given by 

= (l/m) (dP/dn)^ — 5P/3p. Combining these results. 



one immediately obtains the result in Eq. 
bulk viscosity sum rule vanishes there. 
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